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The chiral critical surface is a surface of second order phase transitions bounding the region of 
first order chiral phase transitions for small quark masses in the {m u l i,m s ,jj.} parameter space. 
The potential critical endpoint of the QCD (T, /x)-phase diagram is widely expected to be part of 
this surface. Since for jj. = with physical quark masses QCD is known to exhibit an analytic 
crossover, this expectation requires the region of chiral transitions to expand with fi for a chiral 
critical endpoint to exist. Instead, on coarse N t — 4 lattices, we find the area of chiral transitions 
to shrink with fi, which excludes a chiral critical point for QCD at moderate chemical potentials 
JJ-b < 500 MeV. First results on finer N t = 6 lattices indicate a curvature of the critical surface 
consistent with zero and unchanged conclusions. We also comment on the interplay of phase 
diagrams between the Nf = 2 and Nf = 2 + 1 theories and its consequences for physical QCD. 
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Figure 1: Left: Schematic phase transition behaviour of Nf = 2 + 1 QCD for different choices of quark 
masses (m u j , m s ) at jj. = 0. Right: Measured chiral critical line on an N t = 4 lattice with staggered fermions 




1. Introduction 

The QCD phase diagram has been the subject of intense research over the last ten years. 
Based on asymptotic freedom, one expects at least three different forms of nuclear matter: confined 
hadronic matter (low jx B , T), quark gluon plasma (high T) and colour-superconducting matter (high 
Hb, low T). Whether and where these regions are separated by true phase transitions has to be 
determined by first principle calculations and experiments. Since QCD is strongly coupled on 
scales of nuclear matter, Monte Carlo simulations of lattice QCD are presently the only viable 
approach. 

Unfortunately, the so-called sign problem prohibits straightforward simulations at finite bary- 
on density, thus our expectations for the QCD phase diagram are largely founded on model cal- 
culations. Since 2001, several ways have been designed to circumvent the sign problem in an 
approximate way, all of them valid for \ijT < 1 only [jl|, 0]. Within this range, those methods give 
quantitatively agreeing results for, e.g., the calculation of T c (n) [[§]. Many phenomenologically 
interesting quantities like screening masses, the thermodynamical pressure, quark number suscep- 
tibilities etc. are thus theoretically controlled at moderate quark densities. On the other hand, 
because of the intricate and costly finite size scaling analyses involved, determining the order of 
the QCD phase transition, and hence the existence of a chiral critical point, is a much harder task. 
Here we discuss the problem of determining the order of the finite temperature phase transition by 
lattice simulations in the extended parameter space {m U! d,m s ,T,ix}. 

2. The chiral critical line at ji = 

The order of the QCD finite temperature phase transition as a function of quark masses is 
depicted in Fig. [T], for /J. = (left). In the limits of zero and infinite quark masses (lower left and 
upper right corners), order parameters corresponding to the breaking of the global chiral and centre 
symmetries, respectively, can be defined, and one numerically finds first order phase transitions 
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at small and large quark masses at some finite temperatures T c (m u j,m s ). On the other hand, one 
observes an analytic crossover at intermediate quark masses, with second order boundary lines 
separating these regions. Both lines have been shown to belong to the Z(2) universality class of 
the 3d Ising model 0,^^]. Since the line on the lower left marks the boundary of the quark mass 
region featuring a chiral phase transition, it is referred to as chiral critical line. 

However, the nature of the Nf = 2 chiral transition is far from being settled. Wilson fermions 
appear to see 0(4) scaling [Q], while staggered actions are inconsistent with 0(4) and 0(2) (for the 
discretised theory) [§]. A recent finite size scaling analysis using staggered fermions with unprece- 
dented lattice sizes was performed in [g]. Again, these data appear inconsistent with 0(4)/0(2), 
and the authors conclude a first order transition to be a possibility. A different conclusion was 
reached in [ffO], in which ^QCD was investigated numerically. This is a staggered action modified 
by an irrelevant term such as to allow simulations in the chiral limit. The authors find their data 
compatible with those of an 0(2) spin model on moderate to small volumes, which would indicate 
large finite volume effects in the other simulations. Finally, from universality of chiral models it is 
known that the order of the chiral transition is related to the strength of the E/a(1) anomaly [|i"T|]. 
In a model constructed to have the right symmetry with a tunable anomaly strength, it has recently 
been demonstrated non-perturbatively that both scenarios are possible, with a strong anomaly re- 
quired for the chiral phase transition to be second order Should the chiral transition turn out 
to be first order, the likely modification of Fig. [j] (left) would be the disappearance of the tricritical 
point, with the chiral critical line intersecting the Nf = 2 axis at some finite m u j and being Z(2) all 
the way. 

A convenient observable to locate and identify the second order boundary lines is the Binder 
cumulant 

B 4 (X)^((X-(X)) 4 )/((X-(X)) 2 ) 2 , X = YW- (2-1) 

It has to be evaluated at the (pseudo-)critical coupling j8 c (m,/x), i.e. on the phase boundary defined 
by a vanishing third moment of the fluctuation, ((X — (X)) }\a = 0. In the infinite volume limit, 
B4 — > 1 , 3 for a first order transition or crossover, respectively. At the second order transition, B4 — > 
1.604 dictated by the 3d Ising universality class to which the chiral critical line belongs. On finite 
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Figure 3: Critical surface swept by the chiral critical line as [i is turned on. Depending on the curvature, a 
QCD chiral critical point is present or absent. For heavy quarks the curvature has been determined ^ and 
the first order region shrinks with fi. 



lattices this step function gets smeared out to an analytic function, as shown in Fig. || (left). The 
intersection points of different volumes serve as estimators of the critical value, Z?4(/3 c (jn c ),m c ). A 
scan of B4(j3 c (m),m) in the Nf = 3 theory is shown in Fig. || (right). Also shown is a simultaneous 
fit of all three volumes to a Taylor expansion of the cumulant around the critical point and exploiting 
the theoretically known behaviour under finite size scaling, 

B 4 (m,L) = 1. 604 + bL l/v (m - m c ) + ... (2.2) 



Such a fit allows to check whether the chosen volumes are large enough to be consistent with finite 
size scaling, as well as extracting the universality class from the values of B4 at the intersection 
point and the scaling exponent, v = 0.63 in the case of 3d Ising. 

Following this recipe by fixing m s and then scanning in m Uj d, the chiral critical line has recently 
been mapped out on = 4 lattices [13], Fig. [j] (right). In agreement with expectations, the critical 
line steepens when approaching the chiral limit. Assuming the Nf = 2 chiral transition to be in the 
0(4) universality class implies a tricritical point on the m s -axis, Fig. (left). The data are consistent 
with tricritical scaling [14] of the critical line with m u ^ and we estimate m'" c ~ 2.ST C . However, 
this value is extremely cut-off sensitive and likely smaller in the continuum, cf. Sec.^ 



3. The chiral critical surface 

When a chemical potential is switched on, the chiral critical line sweeps out a surface, as shown 
in Fig. ||. According to standard expectations in the literature [|l4|], for small m u j, the critical line 
should continuously shift with pi to larger quark masses until it passes through the physical point 
at He, corresponding to the endpoint in the QCD (T,pi) phase diagram. This is depicted in Fig. || 
(left), where the critical point is part of the chiral critical surface. However, it is also possible for 
the chiral critical surface to bend towards smaller quark masses, Fig. || (right), in which case there 
would be no chiral critical point or phase transition at moderate densities. For definiteness, let us 
consider three degenerate quarks, represented by the diagonal in the quark mass plane. The critical 
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Figure 4: ji -dependence of the Binder cumulant on the chiral critical line for Nf = 3 (left) [ |15| | and Nf = 
2 + 1 with physical strange quark mass (right), both for N t = 4. 



quark mass corresponding to the boundary point has an expansion 

^ = i + E4-f • (3-D 
m c (0) £ k \itr) 

A strategy to learn about the chiral critical surface is now to tune the quark mass to m c (0) and 
evaluate the leading coefficients of this expansion. In particular, the sign of ci will tell us which of 
the scenarios in Fig. [I] is realised. 

The curvature of the critical surface in lattice units is directly related to the behaviour of the 
Binder cumulant via the chain rule, 

dam r dB<\ ( dBn , C ">) 



d(an) 2 d(a/l) 2 \dam 

While the second factor is sizeable and easy to evaluate, the /I -dependence of the cumulant is 
excessively weak and requires enormous statistics to extract. In order to guard against systematic 
errors, this derivative has been evaluated in two independent ways. One is to fit the corresponding 



Taylor series of B4 in powers of \i /T to data generated at imaginary chemical potential [13, 15], the 



other to compute the derivative directly and without fitting via the finite difference quotient [15], 



Because the required shift in the couplings is very small, it is adequate and safe to use the original 
Monte Carlo ensemble for anf (0),/l = and reweight the results by the standard Ferrenberg- 
Swendsen method. Moreover, by reweighting to imaginary pL the reweighting factors remain real 
positive and close to 1. 

The results of these two procedures based on 20 and 5 million trajectories on 8 3 x 4, respec- 
tively, are shown in Fig. || (left). The error band represents the first coefficient from fits to imaginary 
/i data, while the data points represent the finite difference quotient extrapolated to zero. Both re- 
sults are consistent, and the slope permits and extraction of the subleading pL coefficient, while 
the combination of all data also constrains the sign of the jj, 6 term. After continuum conversion 
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Figure 5: Left: The chiral critical line moves to smaller quark masses with decreasing lattice spacing 
Right: jx 2 -dependence of the Binder cumulant on the chiral critical line for Nf = 3,N t =6 



the result for Nf = 3 is c\ = — 3.3(3), C2 = — 47(20), C3 < [15]. The same behaviour is found for 
non-degenerate quark masses. Tuning the strange quark mass to its physical value, we calculated 
m" ,d (lJ.) with ci = —39(8) and C2 < 0, Fig. f| (right). Hence, on coarse N t = 4 lattices, the region 
of chiral phase transitions shrinks as a real chemical potential is turned on, and there is no chiral 
critical point for /I5 < 500 MeV. Note that one also observes a weakening of the phase transition 



with jit in the heavy quark case [Qj, in recent model studies of the light quark regime Jig , |19j ], as 
well as a weakening of the transition with isospin chemical potential 



4. First steps towards the continuum, N t = 6 

The largest uncertainty in these calculations by far is due to the coarse lattice spacing a ~ 0.3 
fm on N t = 4 lattices. First steps towards the continuum are currently being taken on N t = 6, a ~ 0.2 
fm. At /I = 0, the chiral critical line is found to recede strongly with decreasing lattice spacing 
fll6| , |D| ]: for Nf = 3, on the critical point m K (N t = A)/m K (N t = 6) ~ 1.8. Thus, in the continuum 
the gap between the physical point and the chiral critical line is much wider than on coarse lattices, 
as indicated in Fig. |5] (left). Preliminary results for the curvature of the critical surface, Fig. || 
(right), result in cy = 7(14), — 17(18) for a LO,NLO extrapolation in /I 2 , respectively. Thus the 
sign of the curvature is not yet constrained. But even if positive, its absolute size appears too small 
to make up for the shift of the chiral critical line towards smaller quark masses, and one would 
again conclude for absence of a chiral critical point below jj, B < 500 MeV in this approximation. 
Higher order terms with large coefficients would be needed to change this picture. 

Note that on current lattices cut-off effects appear to be larger than finite density effects, hence 
definite conclusions for continuum physics cannot yet be drawn. A general finding is the steepness 
of the critical surface, making the location of a possible critical endpoint extremely quark mass 
sensitive, and hence difficult to determine accurately. 



5. The interplay between Nf = 2,3 and Nf = 2+1 

Let us finally comment on the importance of understanding both the Nf = 2, 3 as well as their 
connection before concluding anything for physical QCD. Fig. ^ (left) shows the scenario that 
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Figure 6: Two possible scenarios for connected and disconnected triple lines in the m u d = 0-plane. 



is used in the literature when the argument for the expected QCD phase diagram with a critical 
endpoint is made [14]. The tricritical point at pL = 0, which we have discussed here, is expected to 
be analytically connected by a tricritical line, upon varying m s , to the tricritical point at some finite 
/I in the two-flavour case. This picture arises plausibly if the chiral critical surface behaves as in 
Fig. |3| (left). However if, as on our coarse lattices, Fig. || (right) is realised, the situation might well 
be as in Fig. |6| (right). Even if the chiral Nf = 2 theory does feature a tricritical point at finite jU, 
it need not be connected to the chiral critical surface, and nothing follows for the physical point 
without additional information. 

Furthermore, the shrinking of the critical quark masses with diminishing lattice spacing makes 
it likely that a potential tricritical point, Fig. [j], also moves from a large value for the strange quark 
mass, m u j = 0,m' nc ' ~ 2.8J on N t = 4 fll3|], towards smaller values in the continuum limit. In 
particular, it is possible to have a situation for which m u c i = 0,m' s r ' c < m^ hys . In this case the chiral 
critical surface we have been discussing here would not be responsible for a possible critical point, 
regardless of its curvature, but another surface emanating from the 0(4)-chiral limit. Again, both 
of these scenarios depend on conclusively understanding the situation in the Nf = 2 chiral limit. 



6. Conclusions 

The determination of the order of the QCD finite temperature phase transition as a func- 
tion of quark chemical potential is a maximally difficult problem. Besides the sign-problem, the 
strong quark mass and flavour dependence are responsible for potentially rich structures in the 
{ m u,d jj , T, }X } parameter space which are compute-expensive to disentangle due to the required in- 
tricate finite size scaling analyses in the light quark mass regime. The existing results show that we 
need to control various limiting cases of QCD as well as their connection to Nf = 2 + 1 in order 
to understand if there is a critical point in the QCD phase diagram, and to which critical surface it 
belongs. 



Acknowledgements : 

This work is partially supported by Deutsche Forschungsgemeinschaft, project PH 158/3-1. 



7 



Chiral critical surface of QCD 



Owe Philipsen 



References 

[1] O. Philipsen, Eur. Phys. J. ST 152 (2007) 29 [arXiv:0708.1293 [hep-lat]]. PoS LAT2005 (2006) 016 
[arXiv:hep-lat/05 10077]. 

[2] C. Schmidt, PoS LAT2006 (2006) 021 [arXiv:hep-lat/06101 16]. 

[3] P. de Forcrand and S. Kratochvila, PoS LAT2005 (2006) 167 [hep-lat/0509143]. 

[4] F. Karsch, E. Laermann and C. Schmidt, Phys. Lett. B 520 (2001) 41 [arXiv:hep-lat/0 107020]. 

[5] P. de Forcrand and O. Philipsen, Nucl. Phys. B 673 (2003) 170 [arXiv:hep-lat/0307020]. 

[6] S. Kim, Ph. de Forcrand, S. Kratochvila and T. Takaishi, PoS LAT2005, (2006) 166 
[arXiv:hep-lat/05 10069]. 

[7] A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 63 (2001) 034502 [hep-lat/000801 1]; 
Y. Iwasaki et al, Phys. Rev. Lett. 78 (1997) 179 [hep-lat/9609022]. 

[8] C. W. Bernard et al. [MILC Collaboration], Phys. Rev. D 61 (2000) 1 1 1502 [hep-lat/9912018]; 
E. Laermann, Nucl. Phys. Proc. Suppl. 60A (1998) 180. 

[9] M. D'Elia, A. Di Giacomo and C. Pica, Phys. Rev. D 72 (2005) 1 14510 [hep-lat/0503030]. 

[10] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 73 (2006) 074512 [hep-lat/0603021]. 

[11] R. D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984). 

[12] S. Chandrasekharan and A. C. Mehta, [arXiv:0705.0617 [hep-lat]]. 

[13] P. de Forcrand and O. Philipsen, JHEP 0701 (2007) 077 [hep-lat/0607017]. 

[14] A. M. Halasz et al, Phys. Rev. D 58 (1998) 096007 [arXiv:hep-ph/9804290]. 

[15] P. de Forcrand and O. Philipsen, JHEP 0811 (2008) 012 [arXiv:0808.1096 [hep-lat]]. 

[16] P. de Forcrand, S. Kim and O. Philipsen, PoS LAT2007 (2007) 178 [arXiv:071 1.0262 [hep-lat]]. 

[17] G. Endrodi et al, PoS LAT2007 (2007) 182 [arXiv:0710.0998 [hep-lat]]. 

[18] K. Fukushima, Phys. Rev. D 78 (2008) 1 14019 [arXiv:0809.3080 [hep-ph]]. 

[19] E. S. Bowman and J. I. Kapusta, Phys. Rev. C 79 (2009) 015202 [arXiv:08 10.0042 [nucl-th]]. 

[20] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 77 (2008) 1 14503 [arXiv:0712.2625 [hep-lat]]. 



8 



